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Abstract. 

Properties of 2-dimensional generalizations of sine functions that are symmetric or antisymmetric with 
respect to permutation of their two variables are described. It is shown that the functions are orthogonal 
when integrated over a finite region F of the real Euclidean space, and that they are discretely orthogonal 
when summed up over a lattice of any density in F. Decomposability of the products of functions into 
their sums is shown by explicitly decomposing products of all types. The formalism is set up for Fourier- 
like expansions of digital data over 2-dimensional lattices in F. Continuous interpolation of digital data 
is studied. 

1 Centre de recherches mathematiques, Universite de Montreal, C. P. 6128 - Centre ville, Montreal, 
H3C3J7, Quebec, Canada 

2 Department of Physics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical Uni- 
versity, Brehova 7, 115 19 Prague 1, Czech Republic 

3 Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical 
University, Prague, Trojanova 13, 120 00 Prague 2, Czech Republic 

E-mail: jiri.hrivnak@fjfl.cvut.cz, patera@crm.umontreal.ca 



1. Introduction 

The purpose of this paper is to complete and extend [I] by considering the remaining two families of 
special functions and their properties, namely the generalizations of sine functions of two variables, which 
are either symmetric or antisymmetric with respect to permutations of their two variables. They are 
denoted here by sin^ ^^(x, y), where x,y E K and A, \i E N. The functions are of independent interest. 
The paper [lj is devoted to the study of 2-dimensional symmetric and antisymmetric generalizations of 
the common exponential and cosine functions, namely E^^{x,y) and coSp^(#, y). 

In [2], the functions denoted SIN ± and COS ± were introduced for any number of real variables, 
and their properties were studied. Thus, the functions coSp^(#,y) of [lj, as well as the functions 

sin^^(x,y) considered here, are the functions of [2J specialized to 2D. 

Standard trigonometric Fourier decompositions of functions (continuous or discrete) of two variables 
[3] use special functions formed as products of two trigonometric functions, each depending on one 
variable. The variables are measured along two orthogonal axes. The approach undertaken in [lj, and 
extended here, appears to be the only 2D 'trigonometric' alternative to the standard approach in the 
literature. Our expansion functions are also built as products of two trigonometric functions, but the 
two variables are intertwined within each trigonometric function, so that no substitution of variables 
can bring it to the form used in the standard approach. 

Restriction of the functions of [2] to two variables allows us to be more specific about the details of 
their properties, most notable being their discretization and orthogonality, continuous and discrete. In 
particular, analogs of the four types of standard cosine transforms [HO [6], developed in pQ and in here, 
would not be possible without specific description of discrete domains of orthogonality of the expansion 
functions applicable in the four cases. 

Decomposition of all products of pairs of functions into their sums is described here. It was not 
considered elsewhere. There are three types of products of sine functions, sin + sin + , sin + sin - , and 
sin - sin - , with the same arguments x, y. For completeness, we also decompose all products sin^ cos^, 
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product 


sin + sin + 


sin + sin 


sin sin 


sin + cos + 


sin + cos 


sin cos + 


sin cos 


cos + cos + 


cos + COS 


cos COS 


terms 


cos + 


cos - 


cos + 


sin + 


sin - 


sin - 


sin + 


cos + 


cos - 


cos + 



Table 1. Structure of the ten types of products decomposed in the paper. The second 
row shows the functions appearing in all the terms of a decomposition. 



and cos cos . The general structure of decompositions of products is rather interesting. It is summa- 
rized symbolically in Table [TJ 

There are numerous reasons motivating the study of sin^ and cos^ functions in more than one 
variable. One such reason is the ever presence of trigonometric functions in applications of mathematics 
from the elementary to the most sophisticated. Our immediate motivation is to be able to use the 
functions in Fourier-like analysis and interpolation [9j of digital data in 2D. 

We describe four versions of the sine transforms of functions given on lattices. They differ by certain 
small shifts of their arguments which allow orthogonality of the expansion functions to be maintained. 
These transforms correspond to the well known cosine transforms, labeled as /, //, and IV , where 
the expansion functions have separated variables [4J, namely cos(7rmx) • cos(7rm/). The four variants of 
the antisymmetric sine transforms and the four variants of the symmetric sine transforms correspond 
to different boundary conditions at the boundaries of the fundamental domain. Similarly, there are 
different types of ordinary cosine and sine tranforms in ID [6j. 

In Section 2 we consider the sine functions, followed by some general properties of sine and cosine 
functions in Section 3. Discretization of the antisymmetric and symmetric sine transforms are the 
subject of Section 4 together with their interpolations of types I, II, III, and IV. Three remarks are 
found in the last Section. 

2. Continuous sine transforms in R 2 

First we consider the antisymmetric, then symmetric 2D sine transforms. The discretization of 
transforms is described in Section 4. Here the definitions of symmetric and antisymmetric functions 
cos ± (x,y) are recalled because they appear in decompositions of products of sine functions into their 
sums. 

2.1. Antisymmetric sine functions. 

2.1.1. Definitions, symmetries and general properties. 

The antisymmetric sine functions sin^ ^(x,y) : R 2 — > R are defined as follows: 

= sin(7rAx) sin(7r fiy) — sm(7i{ix) sin(7rAy), A,/i,it,j/GR. (1) 

Clearly the functions are continuous and have continuous derivatives of all degrees in R 2 . A few examples 
of functions are shown in Figure [!} Note that instead of the factor 27r, which was used in [2], we use 
the 'half argument tt. 

The following properties of the functions are verified directly from the definition 0: 



sin (A,A)( X '^) = ( 2 ) 

sin (A,o) ( x ' y) = sin (o,/i) ( x ' y) = ( 3 ) 

8in (A,M)^'^ = - 8i V)(^^) ( 4 ) 

sin (-A,,) ( x ' y) = sii v,-M) ( x > y) = - sin (~-A,-M) ( x ' y) = - sin (\^ ( x ' y) ^ 

sin (\,v)( x ' x ) = ( 6 ) 

sin-^(x,0)=sin-^(0,y) = (7) 



sin 



(A,M) 



fay) 



sin(7rAa;) sin(7rAy) 
sin(7r/ix) sin(7r/iy) 



Date: December 1, 2009. 
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sin (A,M)( _x ' y) = sin 7^( x > -v) = ~ sin 7\,A~ x > -y) = - sin 7\.A x > v)- ( 9 ) 



± (A,M) V 



Because of ^ - we consider only sin, A with A > ji > 0. 

In addition, the functions sin^, ^ with k, I € Z have symmetries related to the periodicity of the sine 
function 



and 



(M) v 



sin (fcjJ) (r, s) = 0, r, seZ. 



(10) 



(11) 



The relations Q - (10) imply that it is sufficient to consider the functions sin^ fc ^(x, y) on the closed 
triangle F(Sf) given by its vertices in R : 



(x,y)GF(Sf) = {(0,0), (1,0), (1,1)}, 



(12) 



called the fundamental domain of the affine symmetric group [8]. The relations ^ and (11) imply that 
sm^ kl j(x,y) vanishes on the boundary dF(S^) of the fundamental domain. 

The graphs of a few lowest functions sin^, ^(x, y) in the fundamental domain are plotted in Figure fl 






bm (3,2) bm (4,2) bm (4,3) 

Figure 1. The contour plots of examples of functions sin^, ^(x, y) in the fundamental domain. 



2.1.2. Continuous orthogonality. 

The functions sin^. ^ are pairwise orthogonal on the fundamental domain F(Sif^), 

J ^ ^sm^(x,y)sm^ k,l,k' ,1' e N, k>l,k f >l f (13) 

Any function / : R 2 —> C that is antisymmetric f(x,y) = —f(y,x), periodic f(x + 2r, y + 2s) = 
y), r, 5 E Z, has continuous derivatives and vanishes on dF(S^) can be expanded in the antisym- 
metric sine functions sinT^ A : 

f(x,y)= ^2 c kl^ k ^(x,y), c k i=4 f(x,y)sm~ kl) (x,y)dxdy. (14) 

r Jfe,zeN \ ^H^) 
I fe>i ; 



2.2. Symmetric sine functions. 
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2.2.1. Definitions, symmetries and general properties. 



Two-dimensional symmetric sine functions sin 



form: 



(A,m) ' 



C are defined for A, \i E K. in the following 



sin 



(a,m) 



sin(7rAa:) sin(7rAy) 
sin(7r/ix) sin(7r/iy) 



sin(7rAx) sin(7r/iy) + sin(7r^x) sin(7rAy) 



(15) 



A few examples of functions are shown in Figure [2j Note that instead of the factor 2tt used in [2], we 
use the 'half argument it. 



The following properties of the functions are verified directly from the definition (15): 



sin (A,0) ( X ' 


y) 


= sm+ 4i) (x,y) = 









(16) 




y) 


= sin+ A) 0r,y) 








(17) 


sin (~-A, M )(y> 


x) 


= sin+_ M) (x,y) 


= -sin+ A _ M) (x 


v) = - 


- sin (W x >y) 


(18) 




0) 


= sin (A,M) (0 '^ = 









(19) 




y) 


= sin ^)(y' x )' 








(20) 




y) 


= sin (W x >-y) 




-v) = 


- sin (W x >y)- 


(21) 



Because of (|16|) - (|18|), we consider only such sin^^ with A > \i > 0. 



The functions sin^, ^ with fc,i G Z have symmetries related to the periodicity of sine function: 



sin+ ^ (x + 2r, y + 2s) = sin^ n (x, y) , r, s E Z 



sm 



(r, s) = 0, r, s E Z. 



(22) 
(23) 



The relations (|20|) - (|22|) imply that it is sufficient to consider the functions sin^.^, fc,i G N, fc > / 



on the fundamental domain F(/S| ff ) [S|- 

The graphs of the lowest symmetric sine functions sin^, ^, k, I E {1, . . . , 3}, k > I in the fundamental 
domain are plotted in Figure [2] 






oiii( 22 ) Dm (3,2) Dm (3,3) 

Figure 2. The contour plots of examples of functions smf kl Jx,y) in the fundamental domain. 



2.2.2. Continuous orthogonality. 

The functions sin J, ^ are mutually orthogonal on the fundamental domain F(Slf^), 

I ^ sin (fc,o ( x > I/) sin (fc/,r) ( x ' y) dx d 2/ = ^hk'h' , k.l.k' ,l f ^n,k>l,k f > (24) 
where G^y is defined by 



2 fc = I 

1 otherwise 



(25) 
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Any function / : R 2 — > C that is symmetric f(x, y) — /(y, x), periodic f(x+2r, y+2s) = y), r, s E Z, 
has continuous derivatives and vanishes for x = 1 and y = can be expanded in the symmetric sine 
functions sin^ A : 

f(x,y)= ^2 c ki^ k ^(x,y), ^i^^Gki 1 f(x,y)smf k ^(x,y)dxdy. (26) 

(k,ien\ ^H^) 
I k>i / 



3. Additional properties of 2D trigonometric functions. 

Properties of 2D generalizations of sine and cosine are closely interwoven. Let us first recall definitions 
± ( 

(A,/i) 1 



of the functions cos^ Jx, y) according to [HE]. 



cos± , (x,y)= IZZZi cTZZi = cos(tv\x) cos(tt W ) ± cos(^x) cos(^Ay) , (27) 



3 (A,/i) v 



cos(7rAx) cos(7rAy) 



cos^fix) cos(7T/iy) 



± 



where the upper or lower signs should be taken simultaneously. 

3.1. Laplace and other differential operators. 

An obvious relation between 2D sine and cosine functions arises from the second derivatives, 

d 2 d 2 

■ Sm a a) ( X ' y) = 71-2 A ^ COS t\ a) ( X ' V) ' 7T7T COS t\ a) ( X ' V) = ^ X ^ s[u t\ a) ( X ' V) ' ( 28 ) 



The functions sin^^x, are eigenfunctions of the Laplace operator 



f d 2 . d 2 \ . ± 



as well as of the operator (28) applied twice: 

d 2 d 2 

Functions sin + satisfy the equality 

d . , , x 
7— smTv \(#, x) 0, 

where n is normal to the boundary x — y. 

3.2. Product decompositions. 

Products of two sin^ Jx,y) functions decompose into the sum of cosine functions. Products of 

cos (A y} f unc ^ ons decompose into the sum of cosine functions, while mixed products, sin^ y){ x t v) ' 
cos^, ,Ax,y), decompose into the sum of sine functions, see Table u\ Common trigonometrical identi- 



3 (A',m') v 
ties, 

2sin(pz) sm(qz) = — cos((jp + q)z) + cos((p — q)z), 

2cos(pz) cos(qz) = cos((p + q)z) + cos((p — q)z), (29) 
2sin(^>2:) cos(qz) = sin((p + q)z) + sin((p — q)z), 

are used when decompositions of products of sin^^(x, and cos^^(x,y) are calculated. 

There are altogether ten different types of products to consider. A concise presentation of the ten 
decompositions is found in Table [2j The following examples are intended to illustrate how the actual 
decomposition is obtained from the Table. 

Example 3.1. 

The top left entry of the table identifies the line where the decomposition of products of the form 
sin^^(x,y) • sin^, ^{x,y) is given. For our example, we selected (A,//) = (2,1), (A 7 ,//) = (3,2). The 
second column entry on the same line indicates that all terms of the decomposition are the functions 
cos (q; y) m subsequent eight entries on that line refer to individual terms of the decomposition. 
Each term has a coefficient ±1. The actual sign of a term is shown on the line. The value of a and j3 is 
identified on the two top lines of the column. Thus for the first term of the example we find at the top 
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(A zb u! 


u 

r 


± X') 

_1_ /\ J 










(+,+ 


) (+,-) 






(+,4 


■) (+,-) 


(- 




-,-) 


sin^ 


• sin + 


cos + 


+ 






+ 


+ 








- 


sin + 


• sin - 


cos - 


+ 






+ 




+ 




+ 




sin - 


• sin - 


cos + 


+ 






+ 




+ 




+ 




sin + 


• cos + 


sin + 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


- 


sin + 


• cos - 


sin - 


+ 


+ 


+ 


+ 












sin - 


• cos + 


sin - 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


- 


sin - 


• cos - 


sin + 


+ 


+ 


+ 


+ 












cos + 


• cos + 


cos + 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


- 


cos + 


• cos - 


cos - 


+ 


+ 


+ 


+ 












cos - 


• cos - 


cos + 


+ 


+ 


+ 


+ 













Table 2. Decomposition of all products of functions sin ± (x,y) and cos ± (x,y) into the 
sum of eight such functions. The first column shows the product. The second column 
contains the function appearing in the decomposition terms. Subsequent columns pro- 
vide: (i) the subscripts labeling the decomposition terms (first line); (ii) the pair of signs 
applicable in the subscripts of each term; (iii) the signs in front of each term in the 
decomposition (remaining lines). 



line that a = 2 ± 3 and (3 = 1 ± 2. Which of the two signs is applicable in a and (3 is specified at the 
head of the column of the first term as (+,+), so that cos^^(x,y) = cos^ 3 ^. Subsequent terms are 



identified in a similar way. The symbol (x,y) is omitted to simplify the expression. From (29) follows 
the presence of the coefficient 4 multiplying each product below. 

In this way, some of the subscripts of the terms of the decomposition are negative. A convenient 
convention is to label the terms by non-negative subscripts. Therefore the symmetries of the functions 
need to be used in order to write the subscripts as positive. If a subscript of a term sin^ ^ should be 
0, the sine term vanishes. 

The four specific examples illustrate the decompositions: 

4 sin^j • sin^ 2 j = cos^ - cos^ _^ - cos^^ + cos^ 1 _^ + cos^ - cos^ - cos + cos J _ 2 ) 

= COS (5,3) ~~ COS (5,l) ~~ COS (3,l) + COS (l,l) + COS (4,4) ~~ COS (4,2) ~~ COS (4,0) + COS (2,0) 

4 sin^^ • sin^ 2 ^ = cos^ - cos^ 5 _^ - cos^_ 1>3 j + cos^_ 1 _^ - cos^ + cos^ 4 _ 2 j + cos ( ,4) ~~ cos (o -2) 

= COS (5,3) - COS (5,l) + C ° S r3,l) +° - + COS (4,2) - COS (4,0) + COS (2,0) 

4 sin+ 1} • cos+ 2) = sin+ 3) + sin+ _ 1} + sin+ 13) + sin+ + sin+ 4) + sin+ _ 2) + sin+ 4) + sin+ _ 2) 

= sin J,3) " sin J,i) " sin 5,i) + sin a,D + sin (M) " sin 5,2) +0 + 
4 sin - 25l) • cos+ 2) = sin - ?3) + sin^ 5 _ 1} + sin - ^^ + sin^ _ 1} + sin - ^ + sin^_ 2) + sin^ 4) + sin - 05 _ 2) 

= Sin (5,3) - Sin (5,l) + Sin (3,l) +° + - Sin (4,2) +° + 

4. Discrete transforms 

The four versions of the discrete sine transforms introduced here use a grid of points (x m ,y n ) p], 
extending over finite regions that differ by their boundaries. The grid is defined by the numbers iV, 
T, and b. The number of points of the grid is TV 2 , constant b G [0,1] is the displacement of lattice 
points from their original position, and T determines argument of sine functions, so that the argument 
is equal to y. To make the distinctions between the four versions of the antisymmetric and symmetric 
transforms easy to compare, we first describe the pertinent regions. 
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Define the 'closed square' K[ a ,a'] by 



a,a'\ 



[a, a 7 ] x [a, a 7 ], 



a, a 7 G Z . 



(30) 



For the antisymmetric discrete sine transforms, it is necessary to define a 'partly-open triangle' K^ Q 



K[o,i] :={(x,y)e[0,l}x[0,l}\x>y}. 



!]• 



(31) 



Similarly, for the symmetric sine transforms, we introduce a 'closed triangle' if^y? which contains one 
additional side, described by x = y, in comparison with the triangle for antisymmetric sine transforms: 



K+ A] :={(x,y) £ [0,1] x[0,l]\x>y}. 



(32) 



4.1. Antisymmetric discrete sine transforms. 

The antisymmetric sine functions are closely related to the antisymmetric exponential functions pQ. 
Four types of discrete antisymmetric sine transforms can be derived from the antisymmetric exponential 
transforms. The idea is to suitably extend given functions and then apply the antisymmetric exponential 
transforms from [T]. 

In order to derive discrete antisymmetric sine transforms, we define the following three functional 
extension operators. They extend a complex function defined on the 'partly-opened triangle' ^ to 
functions defined on the 'closed square' K[_l,l] • 

First, let / : K[o,l] - ► C We define its extension E^f : — > C 

(f(x,y) 0<x<L, 0<y<L 

-f(-x,y) -L<x<0, 0<y<L 

-f(x,-y) 0<x<L,-L<y<0 

U(-x,-y) -L < x < 0,-L < y < 0, 

Secondly, let / : i^[o,i] — ► C We define its extension Rf : K[ 0) 2] — > C to the square as follows 

(f(x,y) 0<x<l, 0<y<l 

f(2-x,y) Kx<2, 0<y<l 

f(x,2-y) 0<x<l, Ky<2 

U(2-x,2-y) Kx<2, 1 <y <2. 



E L f(x,y) := < 



(33) 



Rf(x,y) 



(34) 



For the function / : K, 



[o l] ~~ * we define antisymmetric extensions Af : K\o,i\ 

4f : 



C as 



f(x,y) 


-f(y,x) 



x > y 
x = y 
x < y. 



(35) 



In the first two types of discrete antisymmetric transforms, we consider the given functions fa : 
1 j — > C that have zero values for the lines x = 1 and y = 0. In the last two types, we consider 

the given functions fa : ^ — > C that vanish on the line y = 0. The additional conditions on the 
boundary are derived from properties of sine functions on the regions in view. 

The following extensions of a function fa or fa with corresponding values of TV, T and b are substituted 
into the formula (56) from pQ: 



(I) 


EiAfr : 


K l 


-Li] - 


-C, 


where 


N = 


2M, 


T = 


2, 


b = 


1 


(36) 


(II) 


E x Ah : 


K l 


-Li] - 


-c, 


where 


N = 


2M, 


T = 


2, 


b = 


1/2 


(37) 


(III) 


E 2 RAf 2 : 


K l 


-2,2] ~ 


-c, 


where 


N = 


AM, 


T = 


4, 


b = 


1 


(38) 


(IV) 


E 2 RAf 2 : 


K l 


-2,2] ~ 


-c, 


where 


N = 


AM, 


T = 


4, 


b = 


1/2 


(39) 



We notice that E\Afa has zero values on the boundary of if[_i 5 i], on the axes x,y and on the 
lines x = ±y, and E2RAfa vanishes on the boundary of if [-2,2] > on the axes X ^V an d on the lines 
x = ±y. The different conditions of zero values are due to the fact that in the first and second type of 
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antisymmetric sine transforms, we consider sin^^(x,y) in comparison with the third and fourth type, 
where we consider sin^, ^(| , | ). 

Due to (anti)symmetry and zero values of the extended functions on axes x, y, on the lines x = ±y 
and on borders of ^[-1,1] or if [-2,2]? we obtain the final explicit form of the four antisymmetric sine 
transforms: 



AMDST- I. 



M-l . M-l 

4 



(k,l=l\ (m,n=l\ 
I k>l J I m>n J 



where x m = ^, y n = jfc. 

AMDST- II. 

M A1 , M-l 

/II,-/ \ II,- - - f \ II,- 4rf fc,M^/,M j? / \ • - / 

= 2^ C M sm (M)^'^' C M = p 2^ /i(%,ysm W) (%,^ 

(k,l=l\ ^ (m,n=0\ 

\ k>l J I m>n J 

where x m iy n c?m,m 2 and = 1 for k ^ M. 

AMDST- III. 

M-l 

/III,-/ \ in,- • 



M-l 

r fe,z =o\ 
4 M 

= ^2 ^1 ^m,M^n,M/2 (^ra 5 2/n) sm ^ + i ( x m 5 ; 

f m,n=l 1 
\ m>n J 

where x m = ^, y n = ^. 



AMDST- IV. 



M-l M-l 



^M (^,2/) = 5^ C LT' sin ( fc+ i 5 / + i)( X ^)' 4T' = ^2 5^ /2(^m,yn)sin (fc+i ^ +i) (x m ,y n ) 

r/c,z=o\ 2 2 r m,n=0 \ 2 2 7 

wnere x m — ^ , i/ n — ^ 

4.1.1. Example of antisymmetric sine interpolation. 

We proceed as follows: Choose a continuous model function, sample its value on a lattice, develop 
the digital data into the finite series according to one of the transforms above, interpolate the digital 
data, and compare the resulting continuous function with the model function. 

Our model function is the Gaussian distribution 

( x -x') 2 +(y-y') 2 , x 

f(x,v) = e , (40) 

with the following parameters (x , ,y , ) = (0.707, 0.293) and a = 0.079. The function /, restricted to the 

domain F(S^), is shown in Figure pi 

— ... 11 

We calculate the antisymmetric interpolating sine functions of the type AMDST-II, namely ^ 4 ' , 

and ^12 ~ > usm g grids of density given by M = 4, 7, and 12. The interpolating functions, together 



with the sampling grids, are presented in Figure [4j 
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Figure 3. The function / of (40) is plotted over the domain F(Sf £ ). 






M = 4 



M = 7 



M = 12 









Figure 4. Three examples ^ 4 ' - , W~> and ^ 12 ' _ of antisymmetric sine functions of 
the type AMDST-II interpolating the function / in Figure [3j Sampling points for the 
interpolation are shown as small black dots. Figures in the third line depict the difference 

f-*l>N~- 



4.2. Symmetric discrete sine transforms. 

Symmetric sine functions are related to symmetric exponential functions [lj. Four types of discrete 
symmetric sine functions are derived from the four types of symmetric exponential transforms. We 
apply symmetric exponential transforms from pQ to the extensions of a given function. 

In order to derive discrete symmetric sine transforms, we use functional operators £x, i?, and for the 
functions / : — > C, we define symmetric extension Sf : if[o,i] ^ C by the formula 



Sf(x,y) 



f(x,y) x>y 
f(y,x) x<y. 



(41) 



In the first two types of discrete symmetric transforms, we consider the given functions fi : ^ —> C 
that have zero values for the lines x — 1 and y = 0. In the last two types, we consider the given functions 
f2 : K^q 1 j — ► C that vanish on the line y = 0. The additional conditions on the boundary are derived 
from properties of sine functions on the regions in view. 
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The following extension of a function fa or fa with the corresponding values of N, T and b are 
substituted into the formula (65) from [I]: 



(I) 


EiSfx : 


K \ 


-1-1] - 


-c, 


where 


N = 


2M, 


T = 


2, 


b = 


1 


(42) 


(II) 


EiSh : 


K \ 


-1-1] - 


-c, 


where 


N = 


2M, 


T = 


2, 


b = 


1/2 


(43) 


(III) 


E2Sf2 ■ 


K \ 


-2,2] - 


-c, 


where 


N = 


4M, 


T = 


4, 


b = 


1 


(44) 


(IV) 


E2Sf2 ■ 


K \ 


-2,2] - 


-c, 


where 


N = 


4M, 


T = 


4, 


b = 


1/2 


(45) 



We notice that E\Sf\ has zero values on the boundary of ^[-1,1] and on the axes and E2RSfa 
vanishes on the boundary of i^[_2,2] an d on the axes x,y. The different conditions of zero values are 
due to the fact that in the first and second type of symmetric sine transforms, we consider sin^, ^(x, y) 
in comparison with the third and fourth type, where we consider sin^, ^(| , | ). 

Due to (anti) symmetry and zero values of the extended functions on the axes x,y and on borders of 
^[-1,1] or ^[-2,2] j we obtain the final explicit form of the four symmetric sine transforms: 

SMDST-I. 

M-l 

(x, y)= ° l kt sin J,o^' y ^ 

(k,l=l\ 
I k>l J 

4 M_1 

c kl ~ J[{2q ^1 G mn fi (x m , y n ) sinj, ^ (x m , y n ) 

kl r m ,n=l 1 
\ m>n J 

where x m = ^, y n = ^. 

SMDST-II. 

M ^ M-l 

^M + ( x > ?/) = J2 c £ + sin J,/)( x ' ^' c £ + = M 2 G J2 Oi sin+ >0 (x m , y n ) 

{ k>l J I m>n J 



SMDST-III. 



n m+i n+^ 

where x m y n 



M-l 



^m + (^>y) = 5^ 4!/ ,+ sin J + i^ + i)( x ^)' 

(k,l=0\ 2 ' 2 

I fe>Z J 



4 



C I!!' + ~~ M 2 G ^ d m ^Md n ^MGmnh ( x m, 2/rc) sm J, + i ( x m, 2/n) 

^ f m,n=l 1 2 
\ m>n J 



where — m — n 



SMDST-IV. 



m ' y™ ~~ m • 

M-l 



\k+\,i+\y 

(k,l=0\ 

X k>i j 



4 M - x 



4T ,+ - M 2 Ghi G ™/ 2 y™) sm J+i y^) 

™ r m ,n=0 1 
\ m>n J 

n m+i n+ i 

where x m y n -^g 2 -. 

4.2.1. Example of symmetric sine interpolation. 

Consider the Gaussian distribution / defined by (40). We calculate the symmetric sine interpolating 
functions of the type SMDST-II: ^4 I,+ , and ^ 12 ' + - These interpolating functions, together with 

the interpolating grids, are depicted in Figure [5) 
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M = 4 



M = 7 



M = 12 









Figure 5. Three examples , and V>i2 + of symmetric sine functions of the type 
SMDST-II interpolating the function / in Figure |3j Sampling points for the interpolation 
are shown as small black dots. The third line show the difference / — 



5. Concluding remarks 



Decomposition of products of functions described in the article for all pairs of functions sin^ ^ and 
cos (A/i) P aves the way to a wealth of new properties of these functions [10J, such as a large variety 
of trigonometric-like identities, representation of functions as orthogonal polynomials, and recursion 
relations for their construction. 

In the paper, we make no use of a particular property of the four families of functions, which is 
proving useful elsewhere [10J: The functions of each family split into two mutually exclusive congruence 
classes according to the value of the sum of their subscript A + \i mod 2. For example, it can be seen 
in Table [2] that all terms in the decomposition of one product belong to the same congruence class, and 
that the classes add up during multiplication of functions. 

The well known Weyl formula for the character of an irreducible finite dimensional representation of 
semisimple Lie groups is a ratio of two S- functions studied in [7]. The character functions are endowed 
with many properties that are fundamental to the theory of representations in general. The sin^ 
functions in the paper, which depend on two variables, resemble the S- functions of the group SU(3). In 
fact, in the case of one variable, these functions coincide. It is therefore interesting to explore properties 
of functions formed as an analogous ratio of sin^ functions. 
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